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Abstract. 

& ' We investigate infinite-order phase transitions like the Berezinskii-Kosterlitz-Thouless 

. transition observed in a triangular-lattice three-spin interaction model. Based on 

a field theoretical description and the operator-production-expansion technique, we 
perform the renormalization-group analysis, and then clarify properties of marginal 
operators near the phase transition points. The results are utilized to establish criteria 
to determine the transition points and some universal relations among excitation levels 
to characterize the transitions. We verify these predictions via the numerical analysis 
^) , on eigenvalue structures of the transfer matrix. Also, we discuss an enhancement 

of symmetry at the end points of a critical intermediate phase in connection with a 
£SJ . transition observed in the ground state of the bilinear-biquadratic spin-1 chain. 
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Phase transitions and critical phenomena observed in classical spin systems have been 
investigated for a long time. Their theoretical treatments including numerical ones have 
revealed a variety of features, and also have been offering the interfaces to understand 
real materials. At the same time, the universality observed in phase transitions is 
one of the most important concepts. For the two-dimensional (2D) critical systems, 
it is pronouncedly expressed in terms of the conformal symmetry being possessed by 
relevant effective field theories. The central charge c is the widely-known theoretical 
parameter [1]. In the case c < 1, it appears to almost characterize a universality class, 
i.e., a possible set of critical exponents [2]: The discrete systems such as the Ising and 
the three-state Potts ferromagnets show the second-order transitions whose universality 
classes are given by the conformal symmetries with the rational values of c. On the 
other hand, in the case c > 1, there still exist considerable efforts to understand the 
universalities of phase transitions. 




Figure 1. The triangular lattice A and the sublattice dependent numbers (n a , rib, n c ) 
(see the text). The shaded area exhibits the elementary plaquette of A. 



It is widely known that the systems with strong frustrations sometimes exhibit the 
residual entropies and the critical ground states with c > 1. Those of the triangular- 
lattice Ising [31 HI [5], [6] and the square-lattice three-state Potts antiferromagnets [3 [8] 
are the typical ones with c = 1. Further, the Kagome-lattice three-state Potts [H [10] 
and the square-lattice four-state Potts vertex antiferromagnets [TTl [12] were clarified to 
possess the ground states with c = 2 and c = 3, respectively [131 EH]- in addition, the 
dimer, the loop-gas, and the coloring models (some of them can be related to the spin 
models) are other examples to show the criticality with c > 1 [15] . Also, the frustration 
effects can increase the central charge of the finite-temperature criticality as observed in 
the 2D fully frustrated XY model P3] • So they have been gathering the great attentions 
for long time both theoretically and experimentally. 

On another front, multispin interactions appear to include an effect to enhance the 
central charge: The exactly solved Baxter- Wu model consisting of the three Ising-spin 
product interaction is the most basic one [17] ; it shows the second-order transition whose 
universality is the same as that of the four-state Potts ferromagnet [15] . The Ising and 
the four-state Potts criticalities are of c = 1/2 and c = 1, respectively. Therefore, the 
multispin interactions are expected as an another source to bring about larger values of 
c, although they have not been argued frequently in this context. 

In this paper, we investigate the three-spin interaction model (TSIM) introduced 
a long time ago by Alcaraz et al [191 [2D]- Suppose that (k,l,m) denotes three sites 
at the corners of each elementary plaquette (see the shaded area in figure []]) of the 
triangular lattice A consisting of interpenetrating three sublattices A a , Ab and A c , then 
the following reduced Hamiltonian expresses a class of TSIM: 

H = —jfj; Yl C0S (fk + <fl + <fm) ■ (1) 

^ (k,l,m) 

One model parameter, the temperature T, will be measured in units of J/k-Q. The angle 
variables <pk = 27rnfc/p {rik G [0,p — 1]) are located on sites, and define the Z p clock 
variables [2T] (the Baxter- Wu model is contained as its special case of p = 2). As we 
shall take a quick look at them in section [21 the most intriguing ones are the properties 
of an critical intermediate phase theoretically expected for p > 4, and its instabilities 
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to the ordered and the disordered phases. Alcaraz et al derived the vector Coulomb gas 
(CG) representation of the model (20]. Especially, they provided the renormalization- 
group (RG) analysis based on its similarity to the triangular-lattice defect-mediated 
melting phenomenon which is known as the Kosterlitz-Thouless-Halperin-Nelson- Young 
(KTHNY) theory [22l E31 EU ESI [26] . In the previous paper, we also argued the effective 
description of TSIM based on the symmetry properties and the so-called ideal-state 
graph concept by Kondev and Henley [T3l [HI E7]; and then introduced the vector 
dual sine-Gordon Lagrangian density [28]. Since its criticality is of c equal to the 
number of components of the vector field, c = 2 was theoretically expected. We 
performed the numerical calculations to confirm this and further the properties of the 
low-energy excitations. However, a detailed analysis of the model on and around the 
transition points has not been done yet. Here, based on the effective field theory, 
we shall first perform the RG analysis, and derive the RG equations to describe the 
transitions to the ordered and the disordered phases. In both cases, we discuss a mixing 
of marginal operators along a separatrix embedded in the RG-flow diagram because 
the same argument was done for the sine-Gordon model [29], and its importance has 
been widely recognized in the discussions of the Berezinskii-Kosterlitz-Thouless (BKT) 
transitions [22j [30j [31] (for applications to classical systems, see [32], [33] [3U [35]). We 
then clarify excitation spectra characteristic to the phase transitions observed. For 
these purposes, we shall utilize some formulae which require the so-called conformal field 
theory (CFT) data such as the scaling dimensions and the operator-product-expansion 
(OPE) coefficients. Therefore, we provide detailed explanations of the OPE calculations 
among local density operators in our field theory [36J. 

The organization of this paper is as follows. In section EJ according to our previous 
research, we shall explain our Lagrangian density to effectively describe the low-energy 
and the long-distance behaviors of TSIM. The calculations of OPE coefficients necessary 
for the CFT technology, the RG analysis of phase transitions, and the conformal 
perturbation calculations of excitation spectra up to the one-loop order are performed 
there. In section [3] based on the analysis in section El we shall explain our numerical 
calculation procedure (the so-called level-crossing conditions) to determine the transition 
points. We perform numerical diagonalization calculations of the transfer matrix, and 
then provide their estimates. Further, to serve a reliability, we check some universal 
relations among excitation levels observed in finite-size systems — in short, we shall 
perform the level-spectroscopy analysis of TSIM [29]. The section H] is devoted to 
discussions and summary of the present investigation. An enhancement of symmetry 
at the transition points will be pointed out; we shall discuss its connection with a ID 
quantum spin system. For readers' convenience, we shall provide two appendices: In 
appendix A, we summarize the properties of the critical fixed point of our model which 
include the conformal invariance as well as OPE's among basic operators. In appendix 
B, based on appendix A, we provide details in the derivations of some useful relations; 
these will contribute directly to the analysis of the critical phenomena observed in the 
present model. 



Criticality in triangular-lattice three-spin interaction model 



4 



2. Theory 

2.1. Vector dual sine-Gordon model 

Since the symmetry property is the key to understand the criticality and the phase 
transitions, we shall begin with its description. Adding to the translations and the 
space inversions, the model is invariant under the global spin rotations 

p=a,b,c IgA p 

with sublattice dependent integers (see figured]) satisfying the condition n a + rib + n c = 
(mod p) [20]. This symmetry operation — we denote as (n a , nh,n c ) — can be generated 
from two of the following three fundamental operations: 

£ a :(l,p-l,0), R b :(0,l,p-1), and R c : (p - 1, 0, 1). (3) 

Thus, it is referred to as the Z p x Z p symmetry. They satisfy some important relations, 
e.g., RP = i? c i?b-R a = 1- Based on these properties, we introduced the vector dual 
sine-Gordon model in the 2D Euclidean space [28]. Writing the Cartesian components 
of the position vector x in the space as (x,y) (see figured]), then it is defined by the 
Lagrangian density C = Co + C\ + £2 with 

A) = ^£||W*)f, ( 4 ) 

i=x,y 

£1 = 5^ £ :e' M -*«:, (5) 

||M||=pa* 

C 2 = ^ Y :e iN - @ W:. (6) 

||N||=1 

The symbol ": :" denotes the normal ordering and means the subtraction of possible 
contractions of fields between them. We shall employ the same definitions of the fields 
and the vector charges as those in our previous paper: is the dual field to $ and 
is related as iKdiQ = eijdj® (e^ is the antisymmetric tensor). In figure 2 of [28J, we 
explained the so-called repeat lattice 1Z representing the periodicity of $ [13], [H] . Using 
its frame as the Cartesian coordinate, the primitive vectors of 1Z are given by 

ex = (1,0) and e a =(i^). (7) 
Also, the primitive vectors of the reciprocal lattice 1Z* are given by 

The magnetic (electric) vector charge N (M) is quantized in TZ (7Z*) whose contravariant 
(covariant) element is expressed as n a = e a ■ N (m a = e a ■ M), and satisfies the 
condition n a (m a ) G Z. Using these vectors, the periodicities of the fields are given by 
$ = $ + 27rN and © = @ + 27rM. The metric tensors for TZ and TZ* are defined by 
9af3 = e a ■ e /3 and g a/3 = e a ■ e 13 , respectively. They satisfy the condition g^g 1 ^ = 5@ 
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because of the duality relation e a ■ e 13 = 5@. A small calculation using equations d7j) and 
dHJ) yields their explicit forms: 



1 ± 



9ap = ( i I j and ^ = ( | I . (9) 




The squared norm of, say, M is given by ||M|| 2 = m a m a (m a = g al3 mp). In equation 
(jSJ), the lattice constant of 71* is denoted as a* (= ^fg l 

The action So = J d 2 xCo represents an interface model consisting of two kinds of 
massless scalar fields [13l[Tl]. Using the elements of <E», it is rewritten as 

So = J d 2 x |^0 a (x)a,,0 a (x), (10) 

where y/2(j) a = e a ■ <fr {<p a = g a p(j>^) [the element of is also defined by y/29 a = e a ■ 
(9 a = g a ^0p)]. The factor \/2 is for convenience. Then, the two-point function exhibits 
the logarithmic behavior 

(r(x)0' 3 (O)) o = --i 7 ^lnf-) 2 , (11) 



where r and a are the distance between x and on the basal 2D space and the ultraviolet 
(UV) cutoff constant, respectively (• • -) means the average respect to the free part So- 
This implies that the fields themselves cannot represent physical quantities. However, 
as we summarize in appendix A, the current and the vertex operators defined by them 
are the scaling operators, and represent physical quantities. Since the system defined 
by Sq is critical, and possesses the conformal invariance with c = 2, the interface model 
is in a roughening phase if C\ } 2 are both irrelevant. 

The phase locking potential Ci consists of the six terms with the following electric 
vector charges [see figure [2(b)] : 

±pe\ ±pe 2 , and ±p(e 1 + e 2 ), (12) 

whose lengths are all pa*. In the unit cell of TZ, it produces the p 2 potential minima which 
form the triangular lattice as $i oc k = 2nl a e a /p with l a e [0,p — 1] (see figure 1 in [28]). 
and each of which corresponds to one of the p 2 -degenerate states. From the formula, 
equation (1B.3j) . the RG eigenvalue of C\ is given by 2 — 2p 2 /3K on the Gaussian fixed 



point S , so it becomes relevant for K > p 2 /3. Since the Gaussian coupling K stands 
for the stiffness of the interface, it is roughly proportional to the inverse temperature. 
Therefore, C\ can stabilize the flat phase with the long-range order at low temperature. 

Another potential £2 is defined in term of the dual field 0. The vertex operator 
e creates a discontinuity of $ by amount of 27rN around the point x. This topological 

f There arc three 2D spaces: (i) The basal 2D space of A. Since the Cartesian components can be 
used for the position vector x, we employ the alphabetical subscripts, i, j, to specify them, (ii) The 
2D space in which 7Z is embedded and (iii) its dual space in which 1Z* is embedded. Like the case of 
the crystallography, since we have employed the nonorthogonal primitive vectors, equation l[7|). for the 
2D space to which $ belongs, it is necessary to introduce the dual space spanned by equation JHJ. In 
this case, it is convenient to express the vectors as the covariant/contravariant elements; we thus use 
the Greek alphabets, a, /?, 7, as the subscript/superscript, accordingly. 
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Figure 2. The schematic representation of the vector charges in the Lagrangian 
density C: (a) The magnetic vector charges in C 2 [equation (|13[) ] which represent the 
discontinuities of 4?. (b) The electric vector charges in C\ [equation (fT!2|) ] (the p = 6 
case) which bring about the phase locking potential with the p 2 minima. 

defect is necessary to describe the disordered phase at high temperature. According to 
the RG sense, the most relevant terms are sufficient to be included in the potential, 
so the summation is performed only for the following magnetic vector charges with the 
shortest length 1 [see figure [2(a)] : 

±e 1; ±e 2 , and ± (e 1 - e 2 ). (13) 

From the formula equation ( IB. 31) . the RG eigenvalue of £ 2 is given by 2 — K/2 on Sa- 
lt thus becomes relevant for K < 4 and brings about the disordered phase. 

Theoretically, it has been expected that the critical intermediate phase (Tl <T< 
Th) survives for the case p > 4 [20], and also that the point K = (= p 2 /3) [K = Kn 
(= 4)] where £i (£ 2 ) becomes marginal corresponds to T L (T H ) [28]. In [28], we observed 
the existence of the critical intermediate phase for the p = 6 case by the use of the 
numerical method. Further, since the effective field theory possesses the duality nature 
which is not obvious in the lattice model, we checked the validity of our theory on the self- 
dual point deep in the intermediate phase, and also estimated Tl,h semiquantitatively. 
Here, we shall perform the detailed analysis of the system in the vicinity of the two 
phase transition points. 



2.2. Operator product expansions, three-point functions, and remarks 

To utilize the CFT technology, we shall first clarify the relationship among local 
operators which plays an important role in our discussion. Since the intermediate region 
corresponds to the Gaussian fixed line parameterized by K, the so-called A4 operator 

mi eh], 

Kn 2 

^w^-^Eii^wii 2 ( 14 ) 

i=x,y 

which is proportional to Cq and translates the system along the line, is the most 
important one. The two-point function is given by (M(pc)M(0)) = (a/rf so that 
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equation (j!4p defines the truly marginal operator satisfying the normalization condition. 
Adding to this, we define the local operators proportional to £ 1)2 as 

E :e iM '*«:, (15) 

^ 6 ||M||=pa. 

w(x) = ± Yl ■ eiNe(x) : • ( 16 ) 



Nil 



Their two-point functions are (V(x)V(0)) = (a/r) 2xv and (W(x)W(0)) = (a/r) 2xw 
with the dimensions 

2 P 2 A K <rt 

so they are also in the normalized forms. 

First, let us consider the expansion of the operator product V(x)V(0), which 
becomes important for K ~ K^. While there are 36 terms in the double summations 
with respect to the vector charges, say, M and M', the following two cases are enough to 
be taken into account: (i) M + M' = (six terms) and (ii) ||M + M'|| = pa* (12 terms); 
the other 18 terms are irrelevant. After some calculations using the basic relations in 
appendices, we find that the cases (i) and (ii) mainly give M. and V, respectively. We 
then obtain the expression of the OPE as follows: 

v(x)v(o),-^(^) 2 -^(o) + A (? )- v( o ) + .... (18) 

The part "• ■ •" includes the unit operator, the stress tensor as well as less singular terms. 
It should be noted that the second term in the right-hand side (rhs) appears due to the 
triangular-lattice structure of TZ*, which is highly contrasted to the single component 
case, and brings about differences as we will see in the following. The cross-check of 
equation (1181) can be done by performing the another OPE calculation 

A4(x)V(0) ~ (-Y V(0) + ■ ■ • , (19) 



y/2 V 

which exhibits the symmetry property of the OPE coefficients to satisfy, i.e., Cwx — 

Cmvv (— Cvxv)- Now, we can read off the OPE coefficients as follows: 

Xy 2 

Cwm = 7= an d Cvvv = ~~j=- (20) 

v 2 v6 

Next, we shall consider the region near K ~ K^, and derive the OPE of W(x)W(0). 

Using the basic relations in appendices, it proceeds in parallel with the derivations of 

equations (TT5T) and (1T91 . Then, we obtain 

W(x )W (0) c M(0) + -| W(0) + • • , (21) 

A1(x)W(0)^^(^) 2 W(0) + .... (22) 
Thus, the OPE coefficients are given by 

Cwwm — an d Cwww = -y=, (23) 
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where the nonzero Cyvvvvv is again attributed to the triangular lattice structure of the 
repeat lattice 71. 

As one of the consequences of the above OPE calculations, we can fix the three- 
point functions among operators. In contrast to the single component case, we obtain 
the nonvanishing one for the phase locking potentials, e.g, 

7~) ' (24) 

where r^ is the distance between Xj and x^ (the same relation also holds for W). 
This is because three vectors at the angle of 120 degrees to each other [e.g., pe 1 , pe 2 , 
and — ^(e 1 + e 2 ), as visible in figure [2(b)] satisfy the vector charge neutrality condition 
[201 [28] (an extension of the scalar case [371 138]). and this plays an important role in 
the following discussion. 

Lastly, we shall refer to other operators not listed above. The spin degrees of 
freedom is the most basic one, and is defined as Sf. =: e lipk :. In the previous paper, 
based on its response to the spin rotations, equation ([3]), we argued that its sublattice 
dependent expression is given by 

(S a , S b , S c ) = (: e^ 1 ^* :, : e"^'* :, : e"^ * :) , (25) 

whose dimensions are all xs = 2/3K [see equation ( IB. 31) ] [28] . This is one example of a 
general form of the quantities related to the spin degrees of freedom, i.e., 

0( X ;{w M })= ^M:e iM '*: (w M eC). (26) 

||M||=M 

For instance, for S'a(x), the length of vector charges M equals to a*, and the weights 
are given by (ti) e i, w e i +e 2, w e 2, w_ e i, w^ e i^ e 2, w„ e 2 ) = (0,1,0,0,0,0). In the above, 
we observed that the uniform mode whose weights are independent of the direction 
of M is engaged in the Lagrangian density due to its symmetry property [see the 
definition of V in equation (fl5"])]. However, it is also expected that the staggered mode 
with (w pe i,w p ( e i +e 2), iy p e 2 ,^-pei,^-p(ei+e 2 ),w-pe 2 ) = (+1, —1, +1, —1, +1, — 1) plays an 
important role, while we shall not discuss this issue in detail. 

2.3. Renormalization- group equations 

Since the data necessary for the use of the CFT technology have been obtained, we shall 
here perform the RG analysis of our effective field theory. First, we consider the region 
near K ~ K^, where £2 is irrelevant. For convenience, we define the scaling field yo as 

K = (l + y )K L . (27) 

Then, the system can be described by the fixed-point Lagrangian density (or CFT) 
(i.e., the Gaussian part with K^) perturbed by two marginal operators Ai and V as 

£=!£o + £l = 4 + ^ W + ^V(x). (28) 




Figure 3. The schematic RG-flow diagram. The left (right) panel exhibits the flow 
around the multicritical fixed point (2/0,2/1) = (0, 0) [(yo, Up) = (0, 0)] corresponding to 
the transition temperature Th (Tl). The scale of j/o in the left is different from that 
in the right, and the y\ and the y p axes are not on the same plain. The separatrixes 
around the points are given by the dotted lines. 



For a perturbed CFT defined by the Lagrangian density, C gcn , = ^o + X^ -VC/i( x )/27ra 2 , 
where marginal scalar operators O p are normalized as (O (U (x)O i/ (0))q = 5 pv {a/r) A 
((• ■ -)q means the average at the fixed point under consideration), the one-loop RG 
equations are governed by the OPE coefficients: For the change, a — > (1 + d/)a, they are 
given by d\ p /dl = — | ~Y^ V p C^ up \ u \ p (C* up denotes the value on Cq) [39]. In the present 
case, using coefficients of equation ( |20|) at K = K^, we obtain the following equations: 

dy ° (0 - 3*(/) 2 , (29) 



dZ 

-2y (l)y p (l)-y p (l) 2 . (30) 



dl 

Similarly to the BKT RG-flow diagram, these exhibit one separatrix between the ordered 
and the critical phases, i.e., 

y P {l) = -yo{l), (31) 

and one straight flow, y p (l) = 2yo(l)/3, renormalized to the strong-coupling fixed point 
(see the right panel in figure [3]). These are similar to those obtained in the research on 
the triangular-lattice defect melting problem j26j HQ] (see also [20] ). Consequently, we 
can introduce the small parameter t to control the distance from the separatrix as 

y p (l) = -(l + t)y (l) (|t|«l). (32) 

Next, we shall derive the RG equations near K ~ Kn, where C\ is irrelevant. We 
redefine the scaling field yo as 

K = (l+y )K H . (33) 

Then, the system is described by 

£:,£„ + £ 2 = £ H + v^ Mx)+ v^ w(x)i ^ 
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where £q is the fixed-point Lagrangian density for the high-temperature transition. The 
RG equations are similarly obtained as follows: 
dy (l) 



dl 
dyi(0 



-3yi(/) 2 , (35) 
-2y (l)yi(l)-yi(l) 2 . (36) 



dl 

Since these are related to equations (129]) and (!30|) via the replacement (yo,y p ) — > 
(—yo,Vi), one separatrix between the disordered and the critical phases, 

s/i(0 = yo(0, (37) 

and one straight flow, ?/i(Z) = —2y (l)/3, renormalized to the high-temperature fixed 
point, are embedded (see the left panel in figure [3]). Thus, for the same aim, we shall 
introduce the small parameter t as 

yi(0 = (i + *)yo(0 (1*1 <i)- (38) 

Here, we note that Boyanovsky and Holman performed the RG of the vectorial sine- 
Gordon field theory based on the simply-laced Lie algebras pET] . While they provided a 
general argument on properties of operators, and some of them are the same as those 
observed above, we have focused on C corresponding to the specific model, TSIM. 

2-4- Mixing of marginal operators 

According to one of the present authors' discussion for the sine-Gordon field theory, 
linear combinations of marginal operators play an important role [29J. As we see in 
the following, it is true also in the present case. So, we shall consider this issue in 
this subsection. Let us start with the system around the separatrix equation (I3T1) . and 
consider the following quantities: 

AocM+^V and BocV + c 2 M. (39) 

The two real coefficients ci )2 are to be determined from the orthogonality condition 
(^t(x 1 ) J B(x 2 )) = which persists under the renormalization along the separatrix [while 
the normalization conditions, e.g., (^4(xi)^4(x 2 )) = (a/ri 2 ) 4 , are used to determine 
overall constants]. Instead of the correlation function, we consider a more convenient 
quantity [12]: 

F{r 12, a, yo(l),y P (l)) = (^) V(xiMx 2 )>, (40) 

and evaluate F and its response to the change of the cutoff dF/ dl up to the lowest order 
in the coupling constants t/q an d y P - For this, we first expand F with respect to y p as 

■^0*12, a, y Q , y p ) ~ F (r 12 , a, y ) - VQypF^ru, a, y ), (41) 



where 



F =(^) 4 M( Xl )i3(x 2 )) , (42) 
F 1 = (^) 4 / f^|(^( Xl )B(x 2 )V(x 3 )) . (43) 
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We should regularize the UV divergence of the integral over x 3 in equation (I43p by 
excluding two circles of the radius a centered at Xi and x 2 . Explicitly, the integral is 
restricted as 

r H(r 13 -a)H(r 23 -a), (44) 

where H(x) is the Heaviside step function. Noticing 4 — 2xy — 4yo 5 we can rewrite 
equation (1421) as F ~ c\ + c 2 + 4ciy In (f^/a). This exhibits F being almost constant 
F ~ ci + C2, so the condition in the lowest order, 

ci + c 2 = 0, (45) 

should be satisfied. Next, let us consider its response to the change of the cutoff, dF /dl. 
There exist two types of contributions, i.e., (i) a direct one via the cutoff a and (ii) an 
indirect one via the coupling constant y controlled by the RG equations. Since the (3- 
functions ( 1291) and (}30l) only include the second-order terms, we can neglect the latter. 
Then, we obtain 
dF 

^ -4 Cl y . (46) 

The contributions from the change of the coupling constant y p to the response of the 
second term in equation (1411 can be neglected due to the same reason, so we shall 
consider diq/d/ up to the zeroth order in y$. Like the case of the first term, there also 
exist two types of contributions; we can omit the type (ii) contributions. Furthermore, 
as we have already seen in the derivation of equation (j4"6l) . a part of the type (i) 
contributions stemming from the power-of-a factors and giving the 0(?/o) terms can 
be neglected. Consequently, the response is contributed only from the change of a in 
the UV regularization factor 

IT " (^) 4 / M*i) B teM*»)>o J L [H(r 31 - a)H(r 32 - a)] . (47) 

Since the integral is a line one along two circumferences of circles centered at X! and x 2 , 
we can estimate the rhs of equation (|47p by using the asymptotic form of the three-point 
functions [e.g., equation (|24|) ]; the result is the following: 
dFi 

~ -2 [C h MVV (l + Cl c 2 ) + CVvvCi] , (48) 

where C^ vv = —y/2. Consequently, from the lowest-order calculation of the condition 
dF/dl = 0, we obtain the relation ciy + [v3(l + cic 2 ) — ci] y p = 0. On the separatrix 
y p = —y , this is reduced to 

2ci - V3(l + cjca) = 0, (49) 

which, together with equation (1451) . can determine the coefficients. While the quadratic 
equation for ci, c\+2ci/ \/3— 1 = 0, possesses two solutions l/\/3 and —y/3, both of these 
provide an identical description of the operators. Thus, in the following discussion, we 
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choose C\ — \j \/3, and call A and B as the .M-like and the V-like operators, respectively 
(see reference [29]). Their normalized expressions are then given by 




cos t?L sin 
■ sin # L cos t?L 




(50) 



with tan-#L = 1/V3- Here, note the followings: Since the condition to determine the 
mixing angle equation (|49p . is expressed in terms of the OPE coefficients, we can 
recognize it as an appearance of the universal properties of the fixed point Cq. Further, 
in the single component case, the corresponding mixing angle is given by tantf = l/y/2 
[29] . This difference mainly stems from the Cyyv contribution absent in the scalar case. 

Next, let us move on to the region near Th, and consider the system around the 
separatrix equation (137)) . where the following linear combinations are to be determined: 

CocM + diW and V ocW + d 2 M. (51) 

In the same way as the above, the orthogonality condition (C(xi)X>(x2)) = persisting 
under the renormalization along the separatrix determines the real coefficients d\ y2 . Since 
the calculations are performed in parallel with the above case, we can straightforwardly 
derive the equations corresponding to equations fH5|) and (jl9~j) as d\ + d 2 = and 
2di + \/3(l + d\d 2 ) = 0, respectively. In accordance with the above case, the solution 
di = —l/y/3 is chosen, so that C and V are termed as the .M-like and the W-like 
operators, respectively (the difference between A and C is contextually obvious). The 
normalized expressions are then given as follows: 




cos t?H sin ^ H 
■ sin $h cos i?h 




(52) 



with tan^H = — 1 / v^3- Consequently, independently of p, we find a simple relation 
between the mixing angles at the high- and the low-temperature transitions, $l = — 
(=tt/6). 

2. 5. Corrections to finite-size scaling and eigenvalue structures 

We are in position to calculate the renormalized scaling dimensions of operators around 
the fixed points C ' and to discuss the significance of the results above. We shall 
start from the free part defined on an infinitely long cylinder in the x 2 direction with 
a periodicity of L in the X\ direction, and write the partition function using the action 
So, cy i. = J cyl d 2 xC as Z 0)Cy i. = /[d*] e" 5 ^- ex Um T _ 00 'L: e~ T ^. Then, H ,l exhibits 
the Hamiltonian operator associated with the transfer matrix e~ H °' L ; it defines a ID 
quantum system with length L, and is given by Hq,l = L dxi7io(xi) with 



Wo(x) = \ 



71 K 

—TT a (x)iT a (x) H d x (f) a (x)d x <f) a (x) 

K 71 



(53) 



The momentum n a conjugate to the field operator (f> a satisfies [(j) a (x), 7Tp(x')] = i5p5(x — 
x'), and v is the velocity of an elementary excitation. When we writing its eigenvalue 



Criticality in triangular-lattice three-spin interaction model 



13 



and eigenstate as Eq^^ and the CFT provides the finite-size-scaling form of the 
excitation gap as AE 0>L>1/ = EQ Lu —E 0jLjg = 2ttvx u /L, where E 0>L>g and x u are the lowest 
energy and the scaling dimension of the operator corresponding to the state \v) (L is 
supposed to be large enough). Next, we consider the Hamiltonian density corresponding 
to the generic model C gen ,, i.e., H ge nX x ) = ^o( x ) + ^^^{x) /2na 2 . Writing the 
ground-state and the excited-state energies as £x ig and El jU , then we can calculate 
the corrections to scaling within the first-order perturbation as AE^^ = E^ v — E^ g ~ 
^pO^+Su ^vQ]wv) 02]- The parenthesized quantity in the rhs defines the renormalized 
scaling dimension. Using the OPE coefficients, the mixing angle and this formula, 
we obtain the dimensions of the .M-like and the V-like operators near Tl [i.e., near the 
separatrix equation fl32l) ] as 

x A = 2 + 2y (l + h -t\ , (54) 
x B = 2-Gy (1 + 7* ). (55) 



4 r 

Similarly, we obtain those of the A1-like and the W-like operators near Th [equation 
(EE}] as 

xc = 2 - 2y (l + ^tj , (56) 

x v = 2 + 6y (l + ^tj (57) 

(yo was redefined as mentioned above). Since these corrections to scaling are described 
by the OPE coefficients, there are some universal relations among the dimensions. For 
instance, in the present case, we find that 
3x A + x B 

= 2 on y p = -y , (58) 

3x c + x v 

= 2 on yi = y . (59) 

Since the ratio of the level splitting, 1:3, being different from that in the single component 
case [29] is one of features, this can provide a solid evidence of the BKT-like phase 
transition described by the RG equations fl29l and fl30l or equations fl35l) and fl36|) . 



3. Numerical calculations (the p = 6 case) 

In this section, we shall explain our numerical calculations and results to confirm the 
above theoretical predictions. We consider the system on A with M (— > oo) rows of L (a 
multiple of 3) sites wrapped on the cylinder, and define the transfer matrix connecting 
the next-nearest-neighbor rows (see figure [T|). We denote its eigenvalues as X q {L) or 
their logarithms as E q (L) = — |ln|A g (L)| (q specifies a level). Then, the conformal 
invariance in critical systems provides the expressions of the central charge c and the 
scaling dimension the corrections to scaling [HI H3, H6] : 

E g (L) ~ Lf - -^c and AE q (L) ~ ^x q . (60) 




Figure 4. (a) The plot of the level-crossing condition of equation (foTj) (the p = 6 
case). The circles (squares) with the fitting curve exhibit its rhs (lhs). The crossing 
point gives the finite-size estimate of the phase transition point 1/Tl(L) with L = 9. 
The inset shows the extrapolation of the data to the thermodynamic limit, and gives 
1/Tl ~ 1.51. (b) The same plot of equation (|62|) as the panel (a), where the circles 
(squares) with the fitting curve exhibit the rhs (lhs). The crossing point estimates 
l/Ta(L) with L = 9. The inset shows the extrapolation of the finite system data, and 
gives 1/T H ^ 1.05. 



Here, E g (L), AE q (L) [= E q (L) - E g (L)], ( (= 2/y/3), and / correspond to "the ground- 
state energy", "the excitation gap", the geometric factor, and a free energy density, 
respectively. In numerical diagonalization calculations using the Lanczos algorithm, we 
employ two fundamental spin rotations R a and i?b in equation ([3]) as well as the lattice 
translation and the space inversion. This is because the matrix size can be reduced, 
and more importantly discrete symmetries can specify lower-energy excitations. For 
instance, since the spin degrees of freedom on A a transforms as Ra.Sk i— > e l2n ^ p Sk and 
RbSk i— > Sk, the corresponding excitation level can be found in the sector with indexes 
(e l2,r//p ,l), and provides a small scaling dimension xs = 2/3K [28]. Thus, we shall utilize 
also this level for the determinations of the phase transition points. 

First, we consider the system around the separatrix equation (131j) . From equation 
(154]) and the dimension of the sublattice dependent spin, xs — 2(1 — yo)/p 2 , the condition 

p 2 xs = 4 — xa (61) 

is satisfied at t = [i.e., on the separatrix equation (13TT)]. Thus, it can be employed 
as a criterion to determine the low-temperature transition point T^. We perform the 
numerical calculations for the p = 6 case and for the systems up to L = 9. In figure 
11(a), we exhibit temperature dependences of the scaling dimensions, i.e., the both sides 
of equation (IBTj) estimated for the L = 9 site system, and find the level crossing at which 
the condition is satisfied. Therefore, we obtain the finite-size estimate, 1/7l(L), from 
the crossing point. As we show in the inset, the extrapolation of finite-size estimates 
to the thermodynamic limit is performed based on the least-squares-fitting procedure 
1/Tl(L) = 1/Tl + a/L 2 . Then, we obtain the transition point as 1/Tl — 1.51. 
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Figure 5. (a) The check of the universal relations among scaling dimensions at Xl. 
The circles (squares) with the fitting line plot the lhs of equation ([55]) [the difference 
(3p 2 xs — xb)/2] at 1/T = 1.51. The up- and the down-ward triangles show xa and xg, 
respectively, and the diamonds plot p 2 xs (with p = 6). (b) The check of the universal 
relations at Tr. The circles (squares) with the fitting line plot the lhs of equation (|59|) 
[the average (36xs +xp)/4] at 1/T = 1.05. The up- and the down-ward triangles show 
xc and xj), respectively, and the diamonds plot 36x5. 



Second, we consider the determination of Tfi, which can be performed in the parallel 
way to the above. From equation (156|) and the dimension of spin, xs — (1— 2/o)/6, around 
equation (1371) . the condition 

12x 5 = x c (62) 

is found to be satisfied only on the separatrix. In figure IH(b), we provide the same plot as 
figure|4](a), where the circles (squares) with the fitting curve plot the rhs (lhs) of equation 
(|62p . Again, we find the level crossing at which the condition is satisfied. Therefore, we 
can estimate l/Tn(L) from the crossing point. The extrapolation to the limit L — > 00 
is also performed (see the inset), and then we estimate the high-temperature transition 
point as 1/T H ^ 1.05. 

In the previous paper, we roughly estimated the transition points from the behavior 
of the central charge (i.e., deviations from the theoretical value c = 2), and obtained 
V^l,h — 1-5, 1.1, respectively. Thus, our above estimates through the level crossings 
are found to be consistent with the data of the central charge. 

At this stage, it is important to check the universal relations among the scaling 
dimensions. As mentioned, since the relation, equation ( l58i) . must hold between the 
7W-like and the V-like excitations at Tl, we calculate the average (i.e., its lhs) at 
1/T = 1.51. As is shown in figure[5]^a), the estimates converge to the theoretical value 2 
very accurately (see the circles with fitting line), meanwhile the scaling dimensions 
themselves considerably deviate from 2 (see the up- and the down- ward triangles). 
Further, the relation (3p 2 xs — xb)/2 = 2 is also expected to hold at the transition 
point, so we calculate the difference (the p = 6 case), and plot the data in the same 
figure (see the square with the fitting line). Despite the smallness of the system sizes, 
the relation holds within 5% error. These checks can be passed only if the system 
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is at the BKT-like phase transition point, and the numerically utilized levels possess 
the theoretically expected properties. Therefore, these are helpful to demonstrate the 
reliability of our approach and results. 

We perform the same checks for the high-temperature transition. In figure E^b), 
we plot the average [the lhs of equation ( 1591) ] at 1/T = 1.05. We find the excellent 
convergence of the data to 2 in the thermodynamic limit (see the circles with fitting line). 
In addition, another relation (36x5 + %v)/^ = 2 is expected between the dimensions of 
the spin and the W-like operators; we plot the average in the same figure (see the square 
with fitting line). Then, we find the deviation of the limiting value (~ 2.13) from 2. 
This may be due to the following reason: In the thermodynamic limit, the universal 
jump of K(l = oo) from (=4) to occurs at Th, and x$ is inversely depending upon 
K [see equation ( IB. 31 )]. Therefore, xs is sensitive to the temperature. Th was reliably 
estimated from the level crossing, but due to the limitation in the size of systems treated, 
it may include some error which causes the deviation. 

Consequently, we have applied the level-crossing conditions to determine the BKT- 
like transition points, and then we have checked some universal relations among 
excitation levels at the transition points. This strategy (the level spectroscopy) was 
proposed and developed by one of the present authors to analyze the BKT transitions 
in ID quantum spin systems [29]. In that case, the sine-Gordon field theory is relevant to 
the discussion. On the other hand, since the present BKT-like transitions are described 
by the vector sine-Gordon models, we have extended the strategy to be applicable to 
them. Then, we have successfully demonstrated its efficiency through the numerical 
calculations of TSIM. 



4. Discussion and summary 

Up to now, we have concentrated on the properties of the critical intermediate phase: 
It possesses the conformal symmetry with c = 2, and exhibits the transitions to the 
ordered and the disordered phases. While, for the latter, its similarity to the triangular- 
lattice defect melting phenomena was argued in the literature [201 128] , we shall re- visit 
the universality class of the transition, and refer to its relevance to a ground-state phase 
transition observed in a quantum spin chain system. 

In the limit p — > oo, the symmetry of TSIM (Z p x Z p ) becomes the U(l) xU(l), i.e., 
the continuous one, and eliminates the low-temperature ordered phase. The RG-flow 
equations ( 1351) and (1361) still describe the transition to the disordered phase, and enable 
us to analyze the system around the transition point. Based on them, quite recently, 
one of the authors proposed a finite-size-scaling ansatz for the helicity modulus of TSIM 
[32], which has a great relevance to the transition (see also [3S1 US, [50]). The ansatz 
reflects a self- similarity of trajectories with respect to a conserved quantity of the flow, 
and mainly predicts the following: (i) In the disordered phase, the correlation length is 
given by £ oc exp [const /(T— Th)"] with the exponent v = 3/5. (ii) The finite-size-scaling 
function takes a universal value at the transition temperature, which comes from the 
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RG flow along the separatrix y\ = y . While we performed large-scale Monte-Carlo 
(MC) simulations of TSIM in this limit to verify the predictions, here we only mention 
that simulation data exhibit a good agreement with the ansatz (the readers interested in 
the detailed discussions and the calculation results may refer to reference [UJ). Since it 
includes the prediction also in the disordered phase, its confirmation is complementary 
to the present argument of sections [2] and [31 and thus provides another solid evidence 
to support our theoretical description. Simultaneously, the exponent v = 3/5 does not 
agree with the previous research [20], where v = 2/5 was predicted based on the vector 
CG representation and the RG argument on the triangular-lattice defect melting theory 
[25] . Also, in the previous paper [28], the explanation that v takes 2/5 was given based 
on their arguments. But, now we are confident that the exponent should be 3/5, so the 
reason of this discrepancy should be clarified in future. 

Instead of in the 2D classical systems, we can find the same situation in the ground 
state of a ID quantum system. TSIM is invariant under the symmetry group S3 of the 
sublattice permutations, and exhibits the conformal invariance with c = 2 which stems 
from the continuous symmetry in the p — ► 00 limit. Guided by these properties, we are 
led to think of the quantum lattice gas model with three components because it realizes 
the S3 symmetry as permutations of components, and its exact solution shows the c = 2 
criticality [5H (52j [53] . Also, as its extension, the bilinear-biquadratic (BLBQ) spin-1 
chain is widely known, and is defined by the Hamiltonian: 



This model possesses some points where the exact information is available: The Affleck- 
Kennedy-Lieb-Tasaki point [51], the Takhatajan-Babujian (TB) point (6*bt = — 7r /4) 
[551 EE], and the Uimin-Lai-Sutherland (ULS) point (0 ULS = vr/4) [EQ M, [53] . The last 
one which corresponds to the quantum lattice gas model separates the extended critical 
phase (7r/2 > 9 > #uls) [57J [58] and the Haldane phase (#uls > 9 > #tb) [59], and 
it is described by the level-1 SU(3) Wess-Zumino-Witten model. The central charge 
for the former and the correlation length in the latter were calculated as c = 2 and 
£ oc exp [const/ (^uls — #) 3 ^ 5 ], respectively [60], with which the numerical estimations 
agree [571 EE] . According to the analysis by Itoi and Kato, the critical fixed line does 
not exist around the ULS point, so the global RG-flow diagram is considerably different 
from the present one [61]. However, the transition occurs when the system crosses the 
separatrix with the SU(3) symmetry, and if we focusing on the massive region including 
the transition point, the RG feature is seemingly similar to our case. This may be a 
reason why the exponent v takes the same value in both cases. Furthermore, quite 
recently, the SU(iV) self-dual sine-Gordon model consisting of the (N — l)-component 
vectorial fields has been investigated [62],|63]; its relevances to, for instance, the quantum- 
spin chains (including the BLBQ model) and ladders have been indeed clarified [62] . For 
TSIM, we have seen that £ consists of the two current operators j a (z) (a = 1,2). In 
addition, for instance for K = Kn, the potential £2 becomes marginal, and the vector 
charges equation (1T31) in TZ which is isomorphic to the root lattice of the SU(3) Lie 




(63) 



(3,k) 
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algebra provide the six operators vj^(z) with the conformal weight (A, A) = (1,0) (see 
appendix B). Thus, there exist eight chiral current operators, and they may define the 
level-1 SU(3) current algebra [MIE5], as in the case of the Kagome-lattice three-state 
Potts antiferromagnet [131 H3] • From these all, we think that although TSIM is in the 
lower symmetry than that of the BLBQ model, it exhibits a symmetry enhancement at 
the end points of the intermediate phase, and then it may share the same fixed point 
properties with the ULS model while more concrete evidences are desired. 

Lastly, we shall comment on the p = 4 case. Although our theory in section [2] as 
well as the vector CG analysis [20] predicts the intermediate phase for 4 < K < 16/3, 
the previous MC data indicated a sign of the first-order transition between the ordered 
and the disordered phases [19]. However, subsequent studies in computational physics 
revealed that it is in general difficult to distinguish between the weak first-order and the 
second-order transitions just based on MC data [661 IB"?]. Further, there are considerable 
difficulties also in the treatment of the BKT-like phase transitions (e.g., an accurate 
determination of the transition point) by MC methods. Thus, for the p = 4 case, the 
nature of phase transitions does not seem to be established yet. 

As one of the possibilities, other than equations (0J, © and (jHJ), there may exist a 
term inferred from the symmetry consideration, and it might eliminate the intermediate 
phase as speculated before [see figure 7(b) in [20]]. Nevertheless, we think that this issue 
remains as an important future problem. 

To summarize, we have investigated the BKT-like continuous phase transitions 
observed in the triangular-lattice three-spin interaction model (TSIM) based on the 
vector dual sine-Gordon field theory. The basic properties of the local density operators 
(e.g., the scaling dimensions) and their mutual relations (the OPE coefficients) have been 
investigated in detail. Using these CFT data, we have performed the RG analysis of 
phase transitions and the conformal perturbation calculations of the excitation spectra 
up to the one-loop order. Especially, the mixing angles of the marginal operators on the 
separatrixes for the low-temperature and the high-temperature transitions, i.e., t?L,H, 
have been determined and compared to the single component case. Then, we have 
found some universal relations among the renormalized scaling dimensions, which can 
precisely characterize the present phase transitions. Furthermore, we have pointed out 
their importance for the numerical determinations of the phase transition points. To 
check the theory, we performed the numerical diagonalization calculations of the transfer 
matrix of TSIM (the p = 6 case) up to the system size L = 9, and determined the 
transition points as 1/Tl,h = 1.51,1.05, respectively, which was followed by the check 
of the universal relations among the excitation levels. Lastly, we have discussed the 
enhancement of symmetry at the end points of the critical intermediate phase. Based 
on the existence of the eight current operators and the value of the exponent v = 3/5, 
we have argued its relevance to the ground state of the bilinear-biquadratic spin-1 chain. 
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Appendix A. Two-dimensional massless scalars: the operator product 
expansions and the conformal invariance 



The action in equation (flOl) consists of two massless scalars located in the 2D Euclidean 
space. Here, we summarize its basic properties, e.g., the equation of motion, the operator 
product expansions, and the conformal invariance |36j . 

It is convenient to adopt the complex coordinates z, z = x ± iy (the former takes 
the upper sign). When we define (f> a (z, z) = </> a (x), 9 a (z, z) = 6 a (x), d 2 z = 2d 2 x, and 
d, 8 = (d x =F i<9y)/2, then equation (flOl) is expressed as 

S = J d 2 z ^d0 a (z,z)d0 a (z,z). (A.l) 

The classical equation of motion is then 

dd<f) a (z,z) =0, (A.2) 
which exhibits the chiral decomposition of fields, i.e., 

r(z,z), e a (z,z) = ^ir(z)±r(z)}. (a.3) 

In terms of new fields with only holomorphic or antiholomorphic dependence, the action 
is re-expressed as 



So 



j d 2 z l-dMzW a (z), (A.4) 

and their two-point functions are diagonal in the sense that 

(^(^(0)) = -^ln-, (A.5) 

(X 

(r(z)^(0)) = -g^\n-, (A.6) 

a 

which otherwise vanish. These show that ip a and ip a are not the scaling operators. 
However, their derivatives 

j a (z) = iadi) a {z) and ] a (z) = iadip a (z) (A.7) 

exhibit, e.g., (j a {z)j^(0)) = g a/3 (a/z) 2 , so that j a and j a are candidates of those with 
the scaling dimension 1. As usual, this issue can be confirmed by the OPE with the 
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stress tensor which is obtained by the Noether theorem: It is diagonal in the complex 
coordinates, and is given by 

T{z) = ~ : j a (z)j a (z) : and T{z) = \ : j a (z)j a {z) : . (A.8) 

Using the Wick theorem and the Taylor expanding, the OPE's can be obtained as 
follows: 

T(z)j a (0) ~ g)V(0)+ Q W(0), (A.9) 

T(z)] a (0) ~ (|) 2 j a (0)+ (D'a^-CO). (A.10) 

These exhibit that j a and j a are the scaling operators with the conformal weights 
(A, A) = (1,0) and (0,1), respectively. The vertex operators are also important 
examples of the scaling operators; they are introduced by 

v k (z) =. e ifc ^ Q(2) : and v k (z) =: e ^^ (f) : . (A.ll) 

llkll 2 

The two-point functions behave as, e.g., {v k (z)v- k (0)) = (a/z) , where k a is the 

covariant element of a constant vector k and ||k|| 2 = k a k a . Similarly to the above, the 

OPE's of Tv and Tv are given as follows: 

II kll 2 a ^ \ i 

T(z)v k (0) ~ U- Uk (0)+(^J a9uk(0), (A.12) 

r(f)«t(0) ^ H£ (i) 2 ^(°) + (i) 1 ^^ )- (A.13) 

Thus, and ^k are the scaling operators with the weights (||k|| 2 /2,0) and (0, ||k|| 2 /2), 
respectively. 

As we have seen, although the physical quantities possess both the holomorphic 
and the antiholomorphic parts, the OPE's are performed independently in these two 
parts due to the diagonal nature of the two-point functions flA.51) and flA.61) . Therefore, 
for a while, we focus only on the holomorphic part. The OPE of T with itself is given 
by 

T(*)T(0) ~ f (^) 4 + 2 (^) 2 T(0) + QW(0). (A.14) 

Thus, we can read off the central charge as c = 5® = 2, which is equal to the number of 
components of the vector field. The OPE between j a and v k are given by 

j a (z)v k (0)~k a (~yv k (0). (A.15) 

This indicates that j a is the current operator to detect the ath element of the vector 
charge k in the vertex operator. Further, the OPE between two vertex operators plays 
a very important role in our discussion; it can be expressed in the following form: 

v k (z)v_ k ,(0) ~ (^) k ' k : v k „ k ,(0) [l + 0(£)] : . (A.16) 

For the case k ^ k', we can neglect the O (z/a) terms in the rhs. However, for k = k', 
since i'o(z) = 1 by definition, they become important. By expansion, we find 

Z k a j a (0) + \ {k a adj a (0) + ik a3 a (0)] 2 } , (A.17) 



a 



where the O ((z/a) 3 ) terms are dropped. 
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Appendix B. Some useful relations 



In this appendix, we shall derive some useful relations which will be referred to in the 
discussion of section [2j Using equation (LA 71) . the M. operator in equation (TH1) is given 
by 



A*(x) = -—j a (z)f(z). 



(B.1) 



On the other hand, using equation ( 1A.11I) . the vertex operator with the vector charges, 
M and N, is expressed as 

. e i[M.*(x)+N.0(x)] . = Wk ( 2 ) % (^) (B.2) 

with K, K = (K~iM^K + *N)/y/2. From equations (fAl2l and (fATT3|) . we can obtain 
the formula for the scaling dimension of the vertex operator, £m,n = o 

(||K|| 2 +||K|| 2 ). 

It is rewritten as a function of the vectors: 

x m ,n = ^A-^Mf + A||N|| 2 ). (B.3) 



The OPE between the A4 operator and the vertex operator is calculated by using 
equation (IA.15I) as 

K K n 2 



MM 



.[M-*(O)+N-0(O)] . r 



V2 



: c 



i[M-$(O)+N-0(O) 



(B.4) 



where the coefficient is also given by K ■ K = £m,o — #o,n- 

The OPE's between the vertex operators with opposite vector charges are the most 
important part in our calculations. Here, we consider the following quantity: 

Q = I y . eiM.*(x) .. e -iM.#(0) (R5) 

\\M\\=pa* 

where the summation is over the six vectors in equation (fl2l . The product of the 
holomorphic and the antiholomorphic parts gives many terms. Among them, those of 
the first order in the elements of the vector charge M disappear after the summation. 
For the second-order terms, by utilizing the relation, 

> m,„.vn,a = 

6 



m a m p = —g a p, 

||M||=pa* 

we find the following compact expression: 

4p 2 /3K 



(B.6) 



Q 



P 
3K 



z- 2 



-) 2 no) + r-) 



v^A4(0)]}.(B.7) 



Similarly, we can perform the OPE calculation of the following quantity 

n = - V : 

6 ^ 



e iN-e(x) .. e -iN-0(O) 



(B.8) 



||N||=1 



where the summation is over the six vectors given in equation (fT3l) . Like equation ( IB. 61) . 
the relation between the elements of the vector charge N and the metric tensor, 



- E 

6 ^ 



n a n 9 



«/3 



(B.9) 
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is available. So, one can find the expansion 



K 



i 1 + f [© 2t(o) + © 2t(o) + i^r }• ( ri °) 



a 

z 

Consequently, we see that the OPE's include the secondary operators T with (A, A) = 
(2, 0) and T with (0,2) as well as the Ai operator with (1, 1) (see also [68]). In general, a 
rotationally invariant system defined on the plain does not include T and T because they 
possess the conformal spins with length 2. On the other hand, since the A4 operator 
is scalar and proportional to £,q, these results naturally exhibit the renormalizations of 
the Gaussian coupling K caused by the potentials £1,2, and may also indicate that the 
metric tensor has been properly employed to define the fixed-point Lagrangian density 

u 
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